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Fluctuations and Correlations in STAR 



o 
o 

(N 
(N 



> 
(N 
^1" 
O 
(N 
O 

o 

">< 

<D 
i 

=5 



x 



Brijesh K Srivastava, (for the STAR Collaboration) 

Department of Physics 

Purdue University 

West Lafayette, IN 47907, USA 

Received: date / Revised version: date 

Abstract. The study of correlations and fluctuations can provide evidence for the production of the quark- 
gluon plasma (QGP) in relativistic heavy ion collisions. Various theories predict that the production of 
a QGP phase in relativistic heavy ion collisions could produce significant event-by-event correlations and 
fluctuations in, transverse momentum, multiplicity, etc. Some of the recent results using STAR at RHIC 
will be presented along with results from other experiments at RHIC. The focus is on forward-backward 
multiplicity correlations, balance function, charge and transverse momentum fluctuations, and correlations. 

PACS. 25.75.-q Relativistic heavy-ion collisions, 25.75.Gz Particle correlations 



1 Introduction 

The investigation of high energy nucleus-nucleus collisions 
provides a unique tool to study the properties of hot and 
dense matter. The motivation is drawn from lattice QCD 
calculations, which predicts a phase transition from hadro- 
nic matter to a system of deconfined quarks and gluons 
(QGP) at high temperature pQ. The study of event-by- 
event fluctuations provides a novel probe to explore such 
transition in the search for the QGP. It is now possible 
with large acceptance experiments at SPS and RHIC. 



2 Forward- Backward Multiplicity Correlations 

The measurement of particle correlations has been sug- 
gested as a method to search for the existence of a phase 
transition in ultra-relativistic heavy ion collisions [SEE]- 
If the quark-gluon plasma (QGP) is formed in these col- 
lisions, the existence or absence of particle correlations 
could lead to a determination of the presence of partonic 
degrees of freedom. A linear relationship has been found 
in high-energy colliding hadron experiments between the 
multiplicity in a forward r\ region (Nf) and average mul- 
tiplicity in a backward r\ region (iV&)[51|B]: 

(Nb(Nf)) = a + bNf (1) 

The coefficient b is referred to as the correlation coefficient 
and it can be expressed in terms of the expectation value 
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are the backward-forward and 
forward-forward dispersions, respectively. This result 
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Fig. 1. Pictorial representation of Forward-Backward correla- 
tions in pseudorapidity. 



exact and model independent [5] . Short and long-range (in 
rapidity) multiplicity correlations are predicted as a sig- 
nature of string fusion [7|[5]. The short range correlations 
are confined to midrapidity, while long range correlations 
are extended more than 2 units in rapidity. When strings 
fuse, a reduction in the long-range forward-backward cor- 
relation is expected. The existence of long-range multi- 
plicity correlations may indicate the presence of multiple 
partonic inelastic collisions. These correlations arise from 
the superposition of a fluctuating number of strings, such 
that 0: 

(N f N b ) - (N f )(N b ) oc [((« 2 > - (n) 2 )] {N of ){N ob ) (3) 

with (in 2 ) — ( n ) 2 ) the fluctuation in the number of in- 
elastic collisions and (Nof) , (Nob) the average multiplicity 
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Fig. 2. (a) D%f, (b) and (c) b from Au+Au collisions 

at t/snn = 200 GeV, as a function of the STAR reference 
multiplicity, N c h- 

produced from a single inelastic collision. Therefore, 
should be sensitive to the presence of long-range multi- 
plicity correlations. 

We discuss the result on forward-backward multiplic- 
ity correlations from 200 GeV Au+Au collisions for all 
charged particles with px in the range from 0.1 to 1.2 
GeV/c, to ensure a sampling of soft particles only. To 
eliminate short-range correlations, a gap in pseudorapid- 
ity (77) of 1.6 units is considered [4]. The forward pseudo- 
rapidity interval was 0.8 < r\ < 1.0 and the backward was 
— 1.0 < 77 < —0.8 [SJ. 

Fig. 2 shows the results for D^, D'jf, and the correla- 
tion coefficient, b, as a function of N c h for the 8 centrality 
bins. The presence of long-range multiplicity correlations 
are evident from D^j and b. The growth of as a func- 
tion of N c h is consistent with an increasing long-range 
correlation from peripheral to central heavy-ion collisions, 
corresponding to a greater number of fluctuating strings. 

The comparison of and D'jf as calculated from 
the Au+Au data to that from the Parton String Model 
(PSM) [TUlfTTj with the string fusion on or off is presented 
in Figs. 3 and 4, respectively. The PSM (with two string 
fusion) minbias multiplicity distribution closely matches 
that of the corrected Au+Au data. It also describes the 
(pt) enhancement, particle ratios, strangeness production, 
etc., seen in Au+Au collisions at RHIC [TO]. As such, the 
centrality cuts in the PSM correspond to the same per- 
centage of the minimum bias cross section as in the data 



Fig. 3. (a) D% f and (b) D 2 ff as a function of the N par t in 200 
GeV Au+Au collisions, compared to the PSM with no string 
fusion (independent strings). 

and are plotted as a function of the number of participant 
nuclcons in the collision (N part ), calculated from Monte 
Carlo Glauber model [12]. The statistical and systematic 
uncertainties in Figs. 3 and 4 are the same as those de- 
scribed for Fig. 2. Figs. 3 and 4 show good agreement 
in peripheral collisions between data and the independent 
(no string fusion) or collective (with string fusion) model. 
In central Au+Au collisions, D'jj with the independent 
string description deviates from the data, but shows good 
agreement for fusion of two soft strings. This confirms 
the agreement in multiplicity between the PSM (with two 
string fusion) and data. However, there is a large discrep- 
ancy in D^j: for central collisions for both the independent 
and collective PSM, compared to the data. This discrep- 
ancy is greater for the case of independent strings (Fig. 
3). This suggests an additional, dynamical reduction in 
the number of particle sources in central Au+Au colli- 
sions, greater than that provided by the fusion of two soft 
strings (Fig. 4). 

3 Transverse Momentum Distributions and 
String Percolation 

It is postulated that in the collision of two nuclei at high- 
energy, color strings are formed between projectile and 
target partons. These color strings decay into additional 
strings via q — q production, and ultimately hadronize to 
produce the observed hadron yields |13j . In the collision 
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Fig. 4. (a.) Dl f and (b) D 2 ff as a function of the N part m 
200 GeV Au+Au collisions, compared to the PSM with string 
fusion (collective strings). 



process, partons from different nucleons begin to overlap 
and form clusters in transverse space. The color strings 
are of radius ro = 0.20-0.25 fm [13]. The fusion of strings 
to form clusters is an evolution of the Dual Parton Model 
[5] , which utilizes independent strings as particle emitters, 
to the Parton String Model (PSM), which implements in- 
teractions (fusion) between strings |14J . At some point, a 
cluster will form which spans the entire system. This is 
referred to as the maximal cluster and marks the onset 
of the percolation threshold. An overview of percolation 
theory can be found in the following reference [15j . The 
quantity p, the percolation density parameter, can be used 
to describe overall cluster density. It can be expressed as 



P 



S 



(4) 



with N the number of strings, S the total nuclear overlap 
area, and ttTq the transverse disc area. At some critical 
value of p — p c , the percolation threshold is reached. p c is 
referred to as the critical percolation density parameter. 
In two dimensions, for uniform string density, in the con- 
tinuum limit, p c = 1.175 [IB] , 

To calculate the percolation parameter, p, a parameteri- 
zation of pp events at 200 GeV is used to compute the Pt 
distribution 

dN _ a 

dpj ~ (po+PT) n (5) 

where a, po, and n are parameters fit to the data. This pa- 
rameterization can be used for nucleus-nucleus collisions 




Fig. 5. Multiplicity suppression factor, F(p) versus the perco- 
lation density parameter, p. The line is the function F(p) and 
is drawn to guide the eye, not as a fit to the points. The esti- 
mated critical percolation density for 2-D overlapping discs in 
the continuum limit, p c , is shown. 



if one takes into account the percolation of strings by [13] 




(6) 



where Si and S n are the transverse overlap area produced 
by a single and N number of strings respectively. In pp 

collisions at 200 GeV, the quantity (^) = 1.0 ± 0.1, 

due to low string overlap probability in pp collisions. Once 
the pt distribution for nucleus-nucleus collisions is deter- 
mined, the multiplicity damping factor can be defined in 
the thermodynamic limit as [17) 



Hp) 



l-e-P 



P 



(7) 



which accounts for the overlapping of discs, with 1 — e 9 
corresponding to the fractional area covered by discs. 

The percolation density parameter, p, has been deter- 
mined for several collision systems and energies. These 
results have been compared to the predicted value of the 
critical percolation density, p c . If p c is exceeded it is ex- 
pected that the percolation threshold has been reached, in- 
dicating the formation of a maximal cluster that spans the 
system under study. Figure 5 is a plot of the quantity F(p) 
versus the percolation density parameter, (p), for central 
collisions. One can also consider the percolation density 
as a function of centrality in Au+Au collisions. The cen- 
trality expressed in terms of the number of participating 
nucleons (N part ) as found from Monte Carlo Glauber cal- 
culations [5] . More central collisions correspond to greater 
values of N part . Figure 6 shows p as a function of the num- 
ber of participant nucleons in Au+Au collisions at 200 and 
62.4 GeV [TS]. For almost all collision centralities, 200 GeV 
Au+Au exceeds the critical percolation threshold, p c . In 
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Fig. 6. The percolation density parameter, p, as a function 
of collision centrality (N part ) in 62.4 and 200 GeV Au+Au 
collisions. 




0.62 



$-0.6 

< 
V 

0.58 



0.56 



0.54 



0.52- 



Fig. 7. 



"^TaR Preliminary 

-#- Au+Au 

P+P 

-©- d+Au 

Au+Au HIJING 
p+p HIJING 
I .... I 



50 



_L 



100 150 



200 



250 



_L 



300 350 



The balance function widths for p+p, d+Au and 
200 GeV as a function of the 



Au+Au collisions at yj (sjvjv 
number of participating nucleons along with HIJING calcula- 
tions. 



62.4 GeV Au+Au, all except the three most peripheral 
bins exceed p c . 



4 Balance Function 

The balance function is based on the principle that charge 
is locally conserved when particles are produced in pairs 

dung. 



B(A V ) 



(At/) - N + + (Ar,) N- + (At/) - N (Atj)' 



N+ 



N- 



(8) 



where is the number of charged pairs in a given pseu- 

dorapidity range, similarly for iV++. 



iV_+ and N+ 
(N- ) is the number of positive (negative) charged parti- 
cles sum over all the events. Aij = (772 — ?7i| is the width of 
the balance function. If the system exists in a deconfined 
phase for an extended time and the charged pairs are pro- 
duced at hadronization, the pairs will retain more of their 
correlation in rapidity. The balance function may be sen- 
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Fig. 8. v+-,dyn for all charged particles with < 0.5 from 
Au+Au collisions at 20, 130, and 200 GeV and pp at 200 GeV 
as a function of the number of participating nucleons. 



sitive to whether the transition to a hadronic phase was 
delayed, as expected if the quark-gluon phase were to exist 
for a longer time. In Fig. 7 the balance function widths for 
p+p, d+Au and Au+Au collisions at \J (snn) = 200 GeV 
are shown as a function of the number of participants |21j . 



The balance function widths scale with N, 



part 



The widths 



predicted using HIJING calculations are also shown in 
Fig. 7, along with data, and show little dependence on 
centrality and are similar to those measured for pp. The 
measured BZ\(?7) in central Au+Au collisions is consistent 
with trends of model incorporating late hadronization. 



5 Charge Fluctuations 

Combined analysis of fluctuations in, e.g., total and net 
charge for positively and negatively charged particles, as 
well as their ratios can reveal interesting physics. One of 
the quantity related to fluctuations is net charge fluctu- 
ations and is the difference of the number of produced 
positively and negatively charged particles measured in a 
fixed rapidity range, defined as [221123] 




(9) 



where N + and are multiplicities of positive and nega- 
tive particles. The magnitude of the variance, , is de- 
termined both by statistical and dynamical fluctuations, 
V-\ = f-| , s tat + ,dyn- Details of net charge fluctua- 
tions analysis can be found in references [24 . In Fig. 8 the 
v+-,dyn is shown for Au+Au collisions at 20, 130, and 200 
GeV along with the pp data for all charged particles with 
<0.5 and 0.1 < pt < 5 GeV/c as a function of number 

of participants [25) . One can see that ^ yn shows very 

little dependence on incident energy. One must however 
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correct measured ^ yn to account for charge conserva- 
tion |25j . Fig. 9 show the plot of 2/_| ^ yn = _ c iyn + 

4/iVc/j as a function of beam energy. The total charged 
particle multiplicity is given by N c h- Fig. 9 also shows 
the results from other experiments for comparison. The 
net charge fluctuations are smaller than expected based 
on predictions from a resonance gas or a quark gluon gas, 
which undergoes fast hadronization. 



6 p t Fluctuations 

Event-by-event fluctuations of mean p t have been pro- 
posed as a possible signature to search for the phase tran- 
sition [2l)ll27l!28] . Fluctuations involve a statistical compo- 
nent arising from the stochastic nature of particle produc- 
tion, as well as a dynamic component determined by cor- 
relations arising in various particle production processes. 
There are several measures to evaluate mean p t fluctua- 
tions. STAR has used (Ap tii Ap td ) 29\ and Aa pt [30] as 
the measure of the dynamical fluctuations. ^2 pT gives the 
normalized fluctuation in CERES experiment |32J. An- 
other measure, F p t, is defined as a deviation from 1 of 
the ratio of the rms of the event-by-event mean px distri- 
bution in real events to that in mixed events(PHENlX) 
[33 . To characterize transverse momentum correlation, 
the quantity (Apt^Aptj) was calculated. Fig. 10 shows 
(Apt : iAp t: j) for Au+Au collisions at ^/snn = 20, 62, 130 
and 200 GeV. One observes that (Ap t ,iAp t j) decreases 
with centrality. These results are also compared with HI- 
JING calculations as shown in Fig. 10. The values for 
(AptjAptj) predicted by HIJING arc always smaller than 
the data [59] , 

The pseudorapidity and azimuth {q, <fi) bin size depen- 
dence of event-wise mean transverse momentum fluctu- 
ations has been measured in terms of Aa p t [30 , 31j . To 
access underlying dynamics, the corresponding autocorre- 
lations (Ap/ \/Pref) were extracted and are shown in Fig. 
11 for 0-5% central Au+Au collisions. The general form 
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Fig. 11. Autocorrelations, Ap/^/p re f on difference variables 
(jjA, cj>A) for 0-5% central Au+Au collisions. 



of the autocorrelations suggests that the basic correlation 
mechanism is parton fragmentation |31j . 

^2 pT also measures the dynamical fluctuations and 
is proportional to mean covariance of all charged parti- 
cle pairs per event. Fig. 12 shows J2 p t as a function of 
the nucleon-nucleon center of mass energy from SPS to 
RHIC energies. The observed fluctuations at SPS and at 
RHIC are about 1%. F p t is approximately proportional 
to (N) Y, T , where (N) is the mean charged particle mul- 
tiplicity. Fig. 13 shows the magnitude of F pt as a function 
of centrality for Au+Au collisions with p™ ax = 2.0 GeV/c. 
A significant non random fluctuation is seen that appears 
to peak in mid central collisions [33) . 



7 Summary 

Some results from fluctuation and correlation measure- 
ments in STAR have been presented. The forward-backward 
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multiplicity correlation in Au+Au shows collective behav- 
ior and a fusion/percolation approach has been explored 
to understand this. The narrowing of balance function 
width in central collisions is consistent with trends 
predicted by models incorporating delayed hadronization. 
The net charge dynamical fluctuations are found to be 
negative and smaller than the value expected from a reso- 
nance or quark-gluon gas. Transverse momentum correla- 
tions and fluctuations show interesting physics which need 
further exploration. 
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